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Abstract. The time- frequency integrals and the two-dimensional stationary phase method 
are applied to study the electromagnetic waves radiated by moving modulated sources in 
dispersive media. We show that such unified approach leads to explicit expressions for the 
field amplitudes and simple relations for the field eigenfrequencies and the retardation time 
that become the coupled variables. The main features of the technique are illustrated by 
examples of the moving source fields in the plasma and the Cherenkov radiation. It is 
emphasized that the deeper insight to the wave effects in dispersive case already requires 
the explicit formulation of the dispersive material model. As the advanced application 
we have considered the Doppler frequency shift in a complex single-resonant dispersive 
metamaterial (Lorenz) model where in some frequency ranges the negativity of the real part 
of the refraction index can be reached. We have demonstrated that in dispersive case the 
Doppler frequency shift acquires a nonlinear dependence on the modulating frequency of the 
radiated particle. The detailed frequency dependence of such a shift and spectral behavior 
of phase and group velocities (that have the opposite directions) are studied numerically. 

Key words: dispersive media; two-dimensional stationary phase method; electromagnetic 
wave; moving modulated source 
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1 Introduction 

The paper is devoted to applications of time-frequency integrals and the two-dimensional sta- 
tionary phase method for problems of waves propagation from moving sources in dispersive 
media. We consider the electromagnetic fields generated by a moving in a dispersive media 
modulated source of the form 

F(t,x) = a(t)e~ iu ' ot 5(x-x (t)), t € R, x = (x\, X2, X3) G M 3 , 

where coq is an eigenfrequency of the source, a(t) is a slowly varying amplitude, Xo(t) = 
(xoi(t),XQ2{t),xo3(t)) is a vector-function defining a motion of the source, 5 is the standard 
(5-function. 

Some assumptions with respect to sources allow us to introduce a large dimensionless parame- 
ter A > which characterizes simultaneously a slowness of variations of amplitudes and velocities 
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of sources, and large distances between sources and receivers. We obtain a representation of the 
fields as double oscillating integrals depending on the parameter A > 



where F is a complex vector- valued amplitude and S is a real- valued for \u\ large enough phase. 
Generally speaking integral (1.1) is divergent and we consider its regularization which is called 
the oscillatory integral. 

Applying to the integrals (1.1) the stationary phase method we obtain the asymptotics of 
electromagnetic fields for large A > 0. We consider non uniformly moving sources in isotropic dis- 
persive homogeneous media, in a particular case, in the isotropic plasma. Note that this method 
can be apply also for the analysis of waves propagation from moving sources in anisotropic dis- 
persive non homogeneous media and media with negative phase velocity (metamaterials), and 
also for a motion with a velocity larger than a phase velocity of media (the Cherenkov radiation) . 

We would like to note that the asymptotic estimates of one- dimensional integrals are standard 
tools of the electrodynamics (see for instance [17, Chapters 3, 4], [33]) and go back to A. Som- 
merfeld [44], and L. Brillouin [8, Chapter 1]. But in the case of modulated non uniformly moving 
sources the representation of the fields in a form of a one-dimensional integral is impossible. In 
turn, a representation of fields as double time-frequency oscillating integrals with a subsequent 
asymptotic analysis yields effective formulae for both the fields and for the Doppler shifts. In 
particular, it gives new approaches to the Cherenkov radiation (see e.g. [2,3,7,19,21,25] and [29, 
Chapter 14]). In particular, the works [1,2,4-48,50-56] describe properties of the charged par- 
ticle field in different dispersive media including traditional resonant medium [1,2,4-48,51], 
active medium [51], anisotropic medium [5], left-handed medium [14,18,31], and so-called "wire- 
metamaterial" [53]. Some of these works develop a method of analysis of the moving charge field 
using complex function theory methods [5,51]. The papers [13,50] are devoted to investigation 
of the fields of moving oscillators in different media. 

The electromagnetic radiation from moving sources is a classical problem of the electrody- 
namics, and for the isotropic non dispersive media the solution of this problem is given by the 
Lienard-Wiechert potential (see for instance [28, Chapter VIII], [24, Chapter 14]). But the 
Lienard-Wiechert potential is not applicable for dispersive media and our representation is an 
effective tool for the investigation of electromagnetic fields generated by moving sources with a 
variable velocity. 

Note that the above-described method for estimating of the acoustic field generated by moving 
sources in stratified acoustic waveguide has been proposed first in [23], and later on in the 
papers [22,34,35]. 

There is another asymptotic approach to the problems of waves propagation from moving 
sources in dispersive media. It is the ray method in the space of the variables (t, x) (see for 
instance [4,30]). Despite the fact that the ray method is applicable to a wider range of prob- 
lems than the method suggested in the article, its implementation is encountering very serious 
difficulties in solving the ray and transport equations. By contrast the approach proposed here 
leads to simple, having a clear physical meaning, equations for the stationary points, and explicit 
formulae for electromagnetic fields and Doppler shifts. In particular, in the developed approach 
the distinction between the phase and group velocities appear in a natural way. 

The paper is organized as follows. In Section 2 we give auxiliary material concerning the os- 
cillatory integrals and the multidimensional stationary phase method. In Section 3 we consider 
electromagnetic waves propagation from moving modulated sources in dispersive medias. We 
obtain effective asymptotic formulae for the electromagnetic fields, Doppler effects, and retarded 
time. Section 4 devoted to applications obtained in Section 4 formulae. We consider a motion 
with a constant velocity in non dispersive media, electromagnetic field generated by modulated 
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stationary sources in dispersive media, electromagnetic waves propagation from uniformly mov- 
ing sources in the lossless no magnetized plasma (see e.g. [19,24,47,48]). We also formulate the 
equations for the Cherenkov radiation in dispersive medias in terms of representation (1.1) and 
the stationary phase method. 

Further it is emphasized that the deeper comprehension and insight the wave effects in dis- 
persive case already requires the explicit formulation of the dispersive material model. As the 
advanced application of the developed technique in Section 5 it is considered the Doppler fre- 
quency shift in a complex single-resonant dispersive metamaterial (Lorenz) model where in some 
frequency ranges the negativity of the real part of the refraction index can be reached [36, 39, 
42,43,52]. We have demonstrated that in dispersive case the Doppler frequency shift acquires 
a nonlinear dependence on the modulating frequency of the radiated particle. (In dispersive- 
less medium such a function is linear.) The detailed frequency dependence of such shift and 
spectral behavior of phase and group velocities (that have the opposite directions) are studied 
numerically. Discussion and conclusions from our results are found in the last section. 



2 Auxiliary material: stationary phase method 
for the oscillatory integrals 

1°. We use the standard notations for the spaces of differentiable functions: C°°(M n ) is 
the space of all infinitely differentiable functions on M. n , C£°(M. n ) is a subspace of C°°(IR n ) 
consisting of functions bounded with all their partial derivatives, C^(R n ) is a subspace of 
C°°(R n ) consisting of functions with a compact supports. 

2°. We consider integrals of the form 

f f(x)e iS ^dx, (2.1) 

where W 1 B x — > f(x) 6 C m is called the amplitude and the scalar function S is called the 
phase. We suppose that / and S are infinitely differentiable (the existence only of a finite 
number of derivatives is necessary) and satisfy the following conditions: for every multiindex a 
there exists C a > such that 

\d a f{x)\ < C a (x) k , xeR n , (x) = (1 + \x\ 2 y, (2.2) 

(i) S(x) is a real function for |x| is large enough, 
(ii) for every a: \a\ > 2 there exists C a > such that \d a S(x)\ < C a , 
(Hi) there exists C > and p > such that 

|V5(aj)| > C\x\ p 

for | a; | large enough. 

Note that if k > —n integral (2.1) does not exist as absolutely convergent and we need 
a regularization of integral (2.1). Let \ S C^°(W l ), and xi x ) = 1 in a small neighborhood of 
the origin. We set xr{ x ) = x( x /R)- 

Proposition 1 ( [41, Chapter 1]). Let estimate (2.2) and conditions hold. Then there 

exists a limit 

F= lim I X R(x)f( x )e tS(x) dx (2.3) 
independent of the choose of the function x- 
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Proof. We introduce the first-order differential operator L 

Lu(x) = (1 + |V5(a;)| 2 ) _1 (/ - iVS(x) • V)u(as), x € R n . 
One can see that 

£ e iS(x,y) _ e iS{x,y) _ ^.4) 

Let L T be the transpose to L differential operator. Then taking into the account (2.4) and 
applying the integration by parts we obtain 

Fr= I X R(x)f(x)e iS{x) dx = f (Ly( XR (x)f(x))e iS ^dx. (2.5) 

Conditions yield 

\(L' r y(xR(x)f(x))\<C j (x} k -^, 

with the constant Cj > independent of R > 0. Let j > ^jp, then the integral in the right side 
part of (2.5) is absolutely convergent, uniformly with respect to R > 0, and we can go to the 
limit for R — > oo in (2.5). Hence the limit in (2.3) exists, independent of %, and 

F= lim F R = [ ((L T ) j f(x))e iS ^dx, (2.6) 

where j > Mfi. The integrals defined by formula (2.6) are called oscillatory. ■ 

3°. We consider an integral depending on a parameter A > 
f(x)e iXS{x) dx, 



where /, S satisfy condition (2.2), (i)-(m). We say that xq is a non-degenerate stationary point 
of the phase S if 

VS(x ) = 0, 

and 

detS"(aj ) + 0, 

dx dx ) * s * ne Hess matrix of the phase 5 at the point x. 

Proposition 2 (see for instance [6,16]). Let there exist a finite set {xi, . . . , x^} of non- 
degenerate stationary points of the phase S. Then 

N 
3=1 

where 

and sgn S"(xj) is the difference between the number of positive and negative eigenvalues of the 
matrix S"(xj). 
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3 Electromagnetic wave propagation in dispersive media 

The Maxwell equations in dispersive media are obtained by the replacing of the electric and 
magnetic permittivity e, p by the operators e(D t ), p(D t ) where 

-1 /*00 -1 f'OO 

e{D t )u(t) = — / e(oj)u(u;)e iujt du;, n(D t )u(t) = — / p(oj)u{co)e iwt du;, 

to J-oo to J_ 00 

where the Fourier transform 

/>oo 

w(w) = / u(t)e- iu}t dt 



is understood in the sense of distributions. Let 

c(uj) = = 
s/e(u))n{u) 

be a phase velocity, k{uj) = -^-^ be a wave number. 
We suppose (see [29, Chapter IX]) that: 

(i) the functions e(lo),h(uj) are limits in the sense of the distributions of analytic bounded in 
the upper complex half-plane functions; 

(ii) k 2 (co) has a finite number u\ < • • • < ui^ of simple zeros on K, and inf k 2 (oj) > 

for small enough e > 0; 

(Hi) the group velocity v g (co) = yjj^ > for all u G M\[a;i — e, oj^ + e]. 

The system of the Maxwell equations for dispersive medias is 

V xH = e(D t )—+j, 

r) / / 

VxE = p(D t ) — , 

e(D t )V -E = p, 

V -H = 0, 

with the continuity equation 

After the Fourier transform with respect to the time we obtain 

V x H(u) = ie(u)uE{uj) + j{u), (3.1) 

V x E(u) = i/i{u)uH(u), (3.2) 
e(cj)V • jE7(w) = p(w), 

V-H(u) = 0, 

V • j( w ) — iujp(ui) = 0. 

In the standard way system (3.1), (3.2) is reduced to the pair of independent equations 

V x V x E(u) - k 2 (uj)E(uj) = iufi(u)j(u), (3.3) 

V x V x H(u) - k 2 (cj)H(uj) = V x j(u). (3.4) 
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Unique solutions of equations (3.3), (3.4) in domain where k 2 {uj) > is defined by the limiting 
absorption principle. That is 

E(u) = lim E c (u), H(u) = lim HAu), 

where E e (co), H e (uj) are the unique bounded solutions of the equations 

V x V x E e (oj) - (k 2 (uj) + ie)E e (u) = iu fi(u)j (u) , 

V x V x H e [u) - (k 2 (u) + ie)H e (uj) = V x j(u), 

and k(ui) = i^\k 2 (u)\ is a purely imaginary number in domains where k 2 (uj) < 0. For these ui 
equations (3.3), (3.4) have unique decreasing solutions. By using the relations 

V x V x E(u) = -V 2 E(co) + V(V • E(cj)) 

and 

V • E(u) = e~ l {uj)p{u)), 
we reduce (3.3) to the equation 

V 2 E(uj) + k 2 {uj)E{u) = e _1 (w)^(w) - iu(j,{u})j(u) 

= -iuiM(u) (fa) + p|^yV(V • \ = F u . (3.5) 

The similar way we obtain from equation (3.4) that 

V 2 H{uj) + k 2 (u)H(uj) = -V x j(u) = * w (3.6) 

Equations (3.5) and (3.6) are independent and can be used for the definition of E and H. 
Let guj(x) = e 47r |rc| be the fundamental solution of the scalar Helmholtz equation 

&g u {x) + k 2 {Lo) gu] {x) = -8{x), x£R 3 , 

satisfying the limiting absorption principle. Hence the solutions of equations (3.5), (3.6) are 
given as 

1 f'OO 1 f'OG 

E(t,x) = — e^\g u * F u )(x)du>, H(t,x) = — e^\g u * ^)(x)duj. (3.7) 

•J — OO " — oo 

where the convolution is understood in the sense of the distributions 

,ik(u)\x—y\ 



( g *¥)(x) = / — r*(y)dy. 

J R 3 4ir\x - y\ 



Let 

j{t, x) = A(t)v(t)5{x - x (t)), (3.8) 

where 5 is the standard 5-function, v(t) = Xo(t) is a velocity of a source, A(t) is an amplitude 
of the source. Then (3.7) implies that 

1 f / e i(k(u)\x-x (T)\-u(t-r)) 

^iM |x-x,(r) 



H(t,x) = 2-2/ V,, x ( — — ^ v(t) I ^(r)dwdr, (3.9) 
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and 



E(t,x) 



8vr 2 i 



e i(k(u)\x—xo(T)\—u(t—r)) 

A(t)coh(co)(I + k- 2 (co)V x V x •) . --. «(r)da;dr.(3.10) 

\x-Xq(t)\ 



Applying the analyticity of the integrand with respect to co in the upper half-plane C+ we 
deform the line of integration (—00,00) with respect to co into the contour T = (— 00, co\ — e) U 
r' U (co/. + e, +00), where V is located in the upper complex half-plane and bypasses from above 
all singularities {co%, . . . , cok} of the integrand on the real line. 

The phase of the double integrals is 



S(co, t) = k(co) \x — a?o( r )| + ut, 



and 



8S(co,t) \x — Xq(t) 



dco 



+ T, 



where 



V(T,X) 



Vg(u) 

Xq(t) ■ (X - Xq(t)) 
\X - Xq(t)\ 



dS(uj,T) 
dr 



-k{u)v(r, x) + co, 



is the projection of the speed xq(t) on the vector x — xq(t). We suppose that there exists large 
enough R > such that 



inf 

J 2 +T 2 >R 2 



v(t, x) 



Vg(U)) 



1 



>0, 



inf 

J 2 +T 2 >/? 2 



v{t, x) 



c(co) 



1 



> 0. 



Then the phase S in integrals (3.9), (3.10) satisfies the estimate 
|VS(w,r)| > C(\u\ + \r\) 



in the domain { (co, t) G 



: co 



integrals (3.9), (3.10) exist as oscillatory. 



+ t 2 > R 2 } where C = C(R) > for R > large enough. Hence 



3.1 Asymptotic analysis of the fields of moving sources 

Now we introduce a dimensionless parameter A > as 



x 



inf \x 



x (t) 



> 0, 



CO 



where cq is the light speed in the vacuum, > is a characteristic frequency of the problem. 
In what follows we suppose that A is a large parameter, that frequently is a ratio of distance 
between the moving source and the receiver to the field wavelength. Such a distance is much 

more then (j^j f° r an the time. Since a'(t) = jd'(t/X) = O(j) (see equation (3.12)) the 1/A 
in general characterizes the slowness variations of the amplitude a due to slowness variations of 
different parameters of a problem, e.g. charge trajectory, material dispersion, etc. 
In what follows we will suppose that 



A(t) = a(t)e- luJot , 



(3.11) 



where 



a(t) = a(t/X), 



(3.12) 
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a G C£°(M), wo > is an eigenfrequency (a carrier frequency) of the source, 
x (t) = XX (t/X), t G R, 



(3.13) 



A > is a large dimensionless parameter characterizing the slowness of variations of the ampli- 
tude a, and the velocity 

x (t)=X (t/X), teR, 

where the vector-function X (t) G C£°(M) ® C 3 . 

To reduce integrals (3.9), (3.10) to a form containing the large parameter A > we use the 
scale change of variables 

x = XX, t = XT, t = Xl. 

and obtain following representations for magnetic and electric fields 



H X (T,X) 



Ex(Lx) 



1 



8vr 2 A 



rxi 



1 



e i\S{T,X,LU,L) 

\X-X (t)\ 

1 



XW(lu) VxVx ' 



V{l)J doxh, 

\ e iAS(T,X,aj,t,) 



\X-X (L 



(3.14) 

V(i)dudi. (3.15) 



8vr 2 z j r 

The phase S in integrals (3.14), (3.15) at A — > +oo is 

S(T, X, u, l) = k(ui) \X — Xq{l)\ — oj(T — a) — ^o^. 

Note that contributions in the main term of asymptotics of integrals (3.14), (3.15) are given 
by the stationary points of the phase S{T,X,lo, l) located in the domain R\ [oj\ — e,Uk + e}. 
The stationary points of S(T, X , uj, l) with respect to (u, l) for fixed (T, X) are solutions of the 
system 



dS(T,X,uj,i) _ \X -X (l) 

did 

dS(T,X,u,L) 



Vg(u) 



(T — t) = 0, 



di 



-k(u)V(X,T) + (u-uj ) = 0, 



(3.16) 



where 



is the value of the projection of V(l) on the vector X — Xq(l). 

Let oj s = uj s {T,X), i s = l s (T,X) be a non-degenerate stationary point of the phase S. 
It means that (uj s , l s ) is a solution of system (3.16) and 



det S"(T,X,u; s ,l s ) ^ 0, 



where 



S"(T,X,u,l) 



(k"{u)\X-X {i)\ 1-^^ ^ 



\ 



V 9 (U>) 



-k(u 



dV(X,t) 
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is the Hess matrix of the phase S. We denote by sgn5"(T, X,u s , l s ) the difference between the 
number of positive and negative eigenvalues of the matrix S"(T, X,uj s , l s ). 

The contribution of the stationary point (w s , l s ) in the asymptotics of H\(T, X), E\(T, X) 
is given by the formulae (see Proposition 2) 



J\S(T,X,u 3 ,t, s ) 



e 4 



sgn S"(T,X,w a ,i s ) 



Idet S"(T, X, oj s , i s ) 



1/2 



a(i s ) 



x 1 + 



(3.17) 



and 



E\,s(t,x) 



4irXi 



1 



e 4 



sgnS"(T,X,u s ,L a ) 



e i\S(T,X,LO a ,L a ) 

v x v x • — — — V(t« 



Idet S"(T, X,uj s ,l s 



1 1/2 



l + O 



(3.18) 



If the phase S has a finite set of stationary points the main term of the asymptotics of the 
electromagnetic field is the sum of contributions of the every stationary point. 

The expressions (3.17) and (3.18) can be simplified if we are restricted by the terms of the 
order O (j) 

H x ,s(T,X) 

ik(iO s )a(L s ) e ^S(T,X,u la ,L s )+% Sg nS"(T,X M )(y x x | X _ Xq (l s )\ 



4vrA 



and 



E\, a (t,x) 



4ir\i 



|det S"(T,X,u s ,t s )\ 1/2 \X - X (i s )\ 
[co s fi(uj s )V(i s ) - (X7 X V X • V(t s )) \X - X (l.)\] 

e iXS{T,X,ui a ,L a )+ifsgnS"{T,X,ui a , Ls ) 



(3.19) 



\X - Xo(l s )\ \detS"(T,X,u) s ,L s 



,1/2 • 



(3.20) 



Coming back to the variables (t, x) we obtain the following asymptotic formulae 



H s (t,x) ~ — Va, X 
47T 



e iS(t,x,uj a ,T a ) 
\X - Xq(t s 



;V[T, 



a(r s )e 



ifagnS"(t,x,0Js,T s ) 



|det S"(t, x, uj s ,t s 



,1/2' 



(3.21) 



E s (t,x) 



where 



T 



e f sgn S"(t,x,L0 a ,r a ) 

< |detS"'(t )a;j a; s ,r s )| 1/2 ' 

|X-X (T)| 
z - x (f)| = . ■ 



e iS(t,x,ui a ,T a ) 

V X V X - ) -—tv(t s ) 

\x - x {t s )\ 



(3.22) 



A — > oo. 



In formulae (3.21), (3.22) the phase is S(t,x,uj,T) = k(uj)\x — xq(t)\ — u(t — t) — loqt, and the 
stationary points [uj s = ui s (t,x),T s = T s (t,x)) are solutions of the system 



dS(t,x,uj,T) \x — xq(t) 



v 9 (u) 



(t - r) = 0, 
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dS(t, x, uj, t) 
dr 



-k(u)v(x, t) + (uj — ujq) = 0, 



(3.23) 



and 



5 (as, t, uj, r) 



/k*'(u)\x-X (T)\ 

v(x,r) 



Vg[u) 



1 _ v(x,t) \ 



Note that under conditions 



sup 



U(T)I + ^\k"(u)\\x-X (r)\ ) 1. 



sup 

(ui,T)m 2 V w 



dv(x, t) 



Ot 



+ 



Vg(U))\ 



< 1, 



where (To, O) are the scale time and frequency, system (3.23) has an unique solution which can 
be obtained by the method of successive approximations. 

Coming to the variables (t, x) in formulae (3.19), (3.20) we simplify formulae (3.21), (3.22) 



H s (t,x) 



E s {t,x) 



4tt 



|det S"(t,x,uj s ,L s )\ 1/2 \x - x q (t ; 



(3.24) 



a(i>s) 
4ni 



[uj s h(u s )v(t s ) - (V X V X ■ v(t s )) ( | sb - x (t s 



|se-se (t s )| \detS"(t,x,u; s ,T s )\ 1/2 ' 
Example 1. Let xq(t) = (0, vr, H). Then x — Xq(t) = (x%, Xi — vr, x 3 — H), 

\x - x (t)\ = {x\ + (x 2 - vt) 2 + (x 3 - H) 2 ) 1 ' 2 , 
v = (0,^,0). The system for the stationary phase point (ui s (t, x),r s (t, sc)) is 

(x2 + (x 2 -i;t) 2 + (x3-£0 2 ) 1/2 



(3.25) 



-Jb(w) 



u fl (o;) 

W(X2 — I/I") 



(x 2 + (x 2 -t-r)2 + 4) 1/2 



(t - t) = 0, 



+ (w — wo) = 0. 



For applying formulae (3.24), (3.25) we have to use 

<9|se-se (t)| d\x-xo(r) 



(V x xv)\x- x (t) 



dx 3 

X3-H 



,0,v- 



XI 



and 



(V x V x -v)(\x-x (t)\) 



x-x (r)y ' \x-x (t)\ 



X\(X2 — VT) x\ + (x 3 — H) 2 (x 3 — H)(X2 — VT)* 



\x-x (t)\ \x-x (t) 



\X - Xq(t)Y 
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3.2 Doppler effect and retarded time 

Note that for fix point x formulae (3.21), (3.22) can be written of the form 
W{t) = A{t)e iF{ - t \ 

where A(t) is a bounded vector-function, F is a real-valued function such that lim F(t) = oo. 

t— >co 

According to the signal processing theory (see for instance [12]) F(t) is a phase of the wave 
process W(t), and the instantaneous frequency Ui n (t) of the wave process W(t) is defined as 
Win(i) = —F'(t). In our case 

F(t) = S(t,x,uj s (t,x),T s (t,x)) 

= k(u s (t,x)) \x - x (r s (t,x)) | - Lu s (t,x)(t - T s (t,x)) - LU T s (t,x), 

where (uj s (t, x), r s (t, x)) is a stationary point of the phase S. Differentiating of F as a composed 
function we obtain 

dS(t,x,co s (t,x),T s (t,x)) dS(t,x,oj s (t,x),T s (t,x)) du s (t,x) 
-t (t) = 



dt 8uj dt 

dS(t,x,uj s (t,x),T s (t,x)) dr s (t,x) 
~ dr dt ' 

Taking into account that (u s (t,x),T s (t,x)) is the stationary point of S, we obtain that 

It implies that the instantaneous frequency co- m (t) of the wave processes H(t,x), E(t,x) for 
fixed x coincides with oo s (t, x). Hence the instantaneous Doppler effect is 

u) a (t,x) -lo = k(u s (t,x))v(x,T s (t,x)). 

Considering the case k(u) s (t, x)) > we obtain the usual Doppler effect 

v(x, T s (t, x)) > => u s (t,x)>uJo and v(x,T s (t, x)) < ==> co s (t, x) < ujq. 

In the case k(uj s (t,x)) < (metamaterials) we obtain the inverse Doppler effect 

v(x, r s (t, x)) < ==$■ u s (t,x) < coo, and v(x,T s (t,x)) > ==$■ u s (t, x) > ljq. 

It follows from formulae (3.24), (3.25) that T s (t, x) is the excitation time of the signal arriving 
to the receiver located at the point x at the time t. Hence the mode Doppler effect for the time 
(the retarded time) is 

\x — xnirJt, x))\ 

t-T s (t,x)= l - /Vf 'J >0 

because the group velocity v g (uj) > 0. 

4 Applications 

4.1 Moving source in non dispersive medias 

We suppose here that the electric and magnetic permittivity e, fj,, and hence the light speed 
c = —|= are independent of u. We consider a moving source of the form (3.8) where A(t) 



and Xo(t) have form (3.11), (3.13). In this case 

m \ LU\X — ^o( r )| / \ 

b(t, X, U), T) = u{t — T) — UJqT. 
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System (3.23) accepts the form 

\x-x (t)\ _ (t _ r) = Q; _^ u ( a!>r ) + ( w _ wo ) = o (4.1) 
c c 

and the first equation (4.1) independent of to. Note that under subliminal velocity of the source 

sup \v(t)\ < c 

t 

first equation in (4.1) has an unique solution r s = T s (t, x) for every points t, x. Second equation 
in (4.1) implies that 

u> s =u s (t,x) = - v{XtTs) . 

c 

Moreover 

det S"(t, x, co s , r s ) = " (l " V -^f^j 2 , (4.2) 

and 

BgnS"{t,x,u a ,T a ) = 0. (4.3) 

Substituting co s , t s , detS"(u s , r s ), sgn S"(oj s , r s ) from (4.2), (4.3) in formulae (3.21), (3.22) we 
obtain the expression for H(t, x) and E(t, x) 

( e iS(t,ac,uj s ,T s ) . . \ . . 

. V * X { \ X -Mrs)\ V ±s )) 

H(t,x) ^ — 

and 

E(t, x )~ — / sf w 8 > / + —-V X V X - 

47TZ ( l v(x,T s ) } \ k z {0J s ) ) \X — X§[Ts 



for t=\,\x- x (t)\ = A -> oo. 

It should be noted that these formulae are asymptotic simplifications of the Lienard-Wiechert 
potentials [24, Chapter 14]. 

4.2 Modulated stationary source in dispersive media 

Let us consider the electromagnetic field generated by a modulated stationary source of the form 

j(t, x) = A(t)e~ iuj()t 5(x - x )e, u > 0, A(t) = a(t/X), 

where A > is a large parameter, e 6 K 3 is a unit vector, loq > is an eigenfrequency of the 
source. 

Repeating the calculations carried out for obtaining formulae (3.9), (3.10) we obtain 



1 f ( e i{k(w)\x-x \-w{t-T)) 

H(t,x) = —, Va,x : : e A(t)Mt, 

8lT Z J R 2 \ \x-x \ 



I f I e i(k{u))\x-x \-u)(t-r)) \ 

E(t,x) = / A(t)uh(u) (/ + k- 2 (cu)V x V x -) — — — -e dwdr. 

o 71 " % Jm? \ I £C — £Cq | / 
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The further asymptotic analyses of H(t, x), E(t, x) is completely similar to given in Section 4.1. 
The phase S in this case is 

S(t, x, uj, r) = k(u) \x — xq\ — oj(t — r) — ujqt. 

Hence system (3.23) accepts the form 

la; — asol 

T = t ^r— , UJ = UJ , 

Vg{U) 

The phase S(t,x,u},r) has the unique stationary point 

\x-x \ „ fk"(ui Q ) \x - x \ 1\ 

oj s =u , T B = t — - and S (t, x, co , r s ) = . 

v g {u ) V 1 °/ 

Hence det S"(t,x,u)o,T s ) = —1, sgn S(t, x, uo, t s ) = and 

S(t,x,u)o,T s ) = k(u ) \x - x \ - UJ T S . 
It implies that 

I / \x-x \\ f e i(k(u>o)\ X -vo\-uoT s ) \ 
H(t > X) " ^oT J Vx X { I* -sol G ) ' 

1 ' j 4*i V 33 X ) \x-x (r s )\ e 

for the "large" time and distance between the source and the receiver. Note that the retarded 

,. - , fas— ccol 
time is t — J — i — y . 

4.3 Propagation from a moving source in the plasma 

We consider a lossless no magnetized plasma whose the collision frequency equals to zero (see 
for instance [19,47,48]). Hence the constitutive parameters in plasma are 



e(u) = e I 1 - I , /i = Ho, 



Eo, no are the electric and magnetic permittivity of the vacuum, 

2 n Q 2 

co, - 



p meo ' 



where oj p is the plasma frequency, n is the particle density, m, q are the mass and charge of the 
electron. Hence the phase velocity in the plasma is 

c{oj) = 



J- , .9 



the wave-number is 



k(u) = 

c 



UJ 2 — UJ 2 
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and the group velocity is 



v g (uj) = cM/1 



UJ 



UJ 



where c$ is the light speed in the vacuum. 

We consider the electromagnetic field in the plasma generated by a moving source of the 
form (3.8) under conditions (3.11), (3.12), (3.13). The phase S is 



S(t, x, uj, r) 



UJ 2 — uj 2 



CO 



x — xq(t) \ — co(t — r) — u)qt. 



We suppose that ujq > ui p . System (3.23) accepts the form 
(t - r) = 0, 



\X - X Q {T) 



ui 2 — UJ 2 



and under the condition 

K*)l 



CO 



v(x,t) + (uj — UJq) = 0, 



(4.4) 



sup , 



< 1 



system (4.4) has an unique solution (uj s ,t s ) such that co s > uj p . 

The substitution (uj s , r s ) in formulae (3.21), (3.22) gives the expression for the electromagnetic 
field generated by the moving source. 

Example 2. Let v be a constant vector and let v(x,t) = ±|f|. In this case equations (4.4) 
accept the form 



T = t 



\X - Xq(t)\ 



W = Wq± 



UJ 2 — UJ 2 



c VI 



p 
771 



CO 



(4.5) 



where the sign + is taken if the source moving to the receiver and the sign — if the source 
moving from the receiver. We obtain from second equation in (4.5) that 



± 



1 



1 — M 2 



ujq ± MJoj 2 - (1 - M 2 )u\ 



where M = ^ (< 1) is the Mach number. For ui = uif first equation in (4.5) has the unique 
solution Tg. It is easy to see that 



det S"(t,x,uf,T^ 



1± 



coa/1 



(4.6) 



and 



8gaS"(t,x,uf,Tf) = 



(4.7) 
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The substitution (co s ,t s ) and (4.6), (4.7) in formulae (3.21), (3.22) gives the expressions for the 
electromagnetic field 




1 / eiSfafiifri \ a ( r ±) 



1 ± 



M 



u 2 

cqa/1 



for i = |z - x {t)\ = |X ~^ o(T)l , A ^ oo. 
4.4 Cherenkov radiation 

Now we apply the above developed approach to consider a field radiation from a moving with 
a constant velocity v charged particle. Here e is a particle charge, while c{uj) > is a field phase 
velocity in the isotropic dispersive medium. Following [29, Chapter XIV] we use 

p(t, x) = e8(x — vt), j(t, x) = ve5(x — vt). 

We suppose that v = (0,0, v) and x = (xi,X2,xs). For the particle eigenfrequency case coq = 
the phase S reads 

S(t, x, co, r) = k(co) \x — vt\ — co(t — r). 

Hence system (3.23) accepts the form 

dS(t, x.oj.t) co . . , A „. 

V ' ' ' = —r^v(x, t) + co = 0, 4.8 

OT C(C0) 

dS Ml^l = H _ ( t _ T ) = o. (4.9) 

We suppose that the system (4.8), (4.9) has a solution (cj s ,t s ) with a non-trivial frequency 
w s > 0. So, is there exists a pair (co s ,t s ) such that 

v(x, T 8 ) = \v\ cos cp(x, t s ) = c(co s ) > 0, (4-10) 

where cp(x,r) is the angle between v and x — vt. The equation (4.10) can be satisfied if the 
value of the projection of the velocity v on the vector x — vt is positive. 

Let now the pair (co s ,t s ) be an isolated non-degenerated stationary point of the phase S, 

v(x,t s ) = c(co s ), (4.11) 
I x — vr q I 

r s = t- 1 4.12 

and det S"(t, x, co s ,t s ) ^ 0. 

Accordingly to the causality principle the root t s of the equation (4.12) has to be positive. 
It implies the Cherenkov cone condition 

vt-x 3 -\x'\^\(3 2 (co s )-l\ >0, (4.13) 
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where x' = {x\,X2) and 

Vg{0Ja) 

(see, e.g., [1] and references therein). Condition (4.13) connects the time t and position x = 
(x±, X2, X3) of the point where the Cherenkov radiation exists. 

The substitution of (oj s ,t s ) to equations (3.21), (3.22) leads to the following expressions for 
the electromagnetic fields E s (t,x), H s (t,x) (at the point (t,x)) emitted by the moving charge 
with the instantaneous frequency uj s = w s (t,x) > (the Cherenkov radiation) 



H s (t,x 




e -i 



sgn S"(t,x,U! s ,T s ) 



\detS"{t,x,uj s ,r s )\ l/2 

X e(vt-X 3 - \x'\y/\p2(u s ) -1| 



\ ]_ e iS{t,x,tx> 3 ,T 3 ) 

£ 5 (^)~-^(^)(J+^V ffi V,-) v 



e 4 



sgn S"(t,x,u) a ,T a ) 



|det 5"(i, a;, w s , r 
where O(r) is the Heaviside function that is 

0(r) = 



—e^t-xs-l^Vl/? 2 ^)-!!), 




5 Doppler effect in metamaterials. Numerical example 

Further investigations of the field properties in a dispersive medium already requires the know- 
ledge of the spectral properties of the material refraction index. So, for the further we have 
to choose the type of material dispersion, in recent literature normally the Lorenz or Drude 
models, see [14, 31] and references therein. Never the less, it is still difficult to investigate the 
properties of the electromagnetic waves in such dispersive model analytically. (We note that 
recently some semi-analytic methods were developed [53]). As the example of the developed 
approach, below we apply the numerics to study the spectral properties of the Doppler effect 
in a dispersive medium. Here we concentrate on the dispersive metamaterials case where the 
refraction index n can be negative Re(n) < (NIM). (In literature also refer to such a material 
as being left-handed (LH) material). We consider such a medium characterized by a (relative) 
permittivity e(u>) and a (relative) permeability h(uj), both of which are complex functions of 
frequency u and the refraction index n(u) satisfying the relations [15,57] 

n(u) = vWMe i[ ^ (w)+ ^ (u)l/2 . (5-1) 

In order to allow a frequency dependence of the refractive index n, let us restrict our attention 
to a single-resonance permittivity 

e(u) = 1 + 2 i e , (5.2) 

and a single-resonance permeability 



2 



/i(w) = l + -a , (5.3) 



4m - - 100J„ 
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Figure 1. (Color online.) Frequency dependence (a) metaterial refractive index n{uj) in the frequency 
interval with Ren(w) < 0; (b) phase velocity v p (u>); (c) group velocity v g (u>), where w = 2nf. 

where wp e , wp m are the coupling strengths, WT e > ^Tm are the transverse resonance frequencies, 
and 7 e , 7 m are the absorption parameters. Both the permittivity and the permeability satisfy the 
Kramers-Kronig relations [15,57]. The following typical parameters of metamaterial were used 
for our numerics: f Pe = 298.42 THz, 7e = 0.04 THz , / Te = 409.82 THz, f Pm = 171.09 THz, 
7 m = 0.04 THz, / Tm = 397.89 THz. 

Fig. 1(a) shows the frequency dependence of NIM refractive index n(u) (u = 2irf). Here 
the permittivity e{oj) and the permeability /i(uj) being respectively given by equations (5.2) 
and (5.3) in the frequency interval from 410 THz to 432 THz where Ren(w) < 0. It is worth 
noting that the negative real part of the refractive index is typically observed together with 
strong dispersion, so that absorption cannot be disregarded in general. However, in a recent 
experiment [56], it was demonstrated that the incorporation of gain material in a metamaterial 
makes it possible to fabricate an extremely low-loss and active optical devices. Thus, the original 
loss-limited negative refractive index can be drastically improved with loss compensation in 
the visible wavelength range. Following this result we will neglect the imaginary part Im(n). 
Fig. 1(b) shows corresponding frequency dependence of the phase velocity v p (tu) that is negative 
in considered frequency range, and Fig. 1(c) exhibits the frequency dependence of the group 
velocity v g {uo). Further for numerics we renormalized the variables for r and x with the scales Iq/c 
and Iq respectively, velocities are normalized with the vacuum light velocity c, and Iq = 75 nm 
is the typical spatial scale used at the metamaterial experiments [56]. To seek for simplicity, 
further we consider with details a case when the position of observer x and the trajectory of 
a source v are in the same plane (X3 = and H = 0). In this case the geometry becomes 2D 
one and the solution to equations (3.23) for uj reads 
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where p 2 = n(ui) 2 v 2 (x 2 — vt) 2 , and r = y x 2 + [x 2 — vt) 2 . The solution of equations (3.23) 
for r can be written as 



x 2 v - v 2 t - Vgi + q 2 

-,2_„2 ' ( b - b > 



v ~ V 9 



where qi = (t> 2 — v 2 ) x\ , and q 2 = v 2 {x 2 — vt) 2 . In simplest ID situation with x\ = 
from (5.4), (5.5) we have two equations 

w x 2 -v g t 

w = -— — — , t= y -. 5.6 

1 ± n(oj)v v — v g 

First formula in equation (5.6) describes the well-known Doppler effect, when the shift of ui de- 
pends on the nv signum. For conventional materials with n > (in this case we have to choose 
upper + sign) the received frequency uj is lower (compared to the emitted frequency uiq) during 
the source approach v > 0, and u) is higher during the source recession v < 0. However the situa- 
tion becomes more complicated in dispersive metamaterials with negative refraction index (NIM) 
where Re(n(w)) < with the refraction index depending on the frequency n = n(uj) were both 
phase and group velocities are the frequency functions. In dispersive medium the first relation 
in (5.4) and (5.5) become coupled equations. Even in a simple ID geometry the first equation for 
frequency u in (5.6) already requires the spectral details of the medium refraction index n(uo). 
Moreover, in the second equation in (5.6) the group velocity v g becomes the frequency function, 
thus, the retardation time r will be different for different source eigenfrequency coq. 

In our numerical simulations we search for the solution of equations (5.1)-(5.6) for u and r at 
the fixed position x±, x 2 and the time t of an observer. First we studied more simple ID dispersive 
NIM, equation (5.6). This equation is solved by the standard numerical methods [38,49], and 
resulting dependencies ares shown in Fig. 2. Fig. 2(a) shows the shifted frequency u (cj = 2nf) 
vs the source frequency coq for area where n(u) < 0. Further such a calculated dependence 
uj = u}(uq) allows us to evaluate the group velocity v g {oS) and, finely the retardation time 
r = t(u)q) in Fig. 2(b). We observe from Fig. 2(a), (b) that both dependencies u) = lo(uiq) and 
r = t(wo) have pronounced nonlinear shape that certainly is determined by a dispersive spectrum 
of used metamaterial. Let us remind that in dispersiveless case both Dopple's frequency shift 
and the retardation time depend on the particle velocity v only. 

The 2D geometry is more complicated. In this situation we have to solve numerically the 
equations (5.4), (5.5) (with added material relations (5.1)— (5.3)) that becomes a strong nonlinear 
system. It is worth to note that in this case the Newton's method for solving nonlinear equations 
has an unfortunate tendency to wander off [38, 49] if the initial guess is not sufficiently close 
to the root. In order to evaluate the solution of such a system the globally convergent multi- 
dimensional Newton's method was applied [38,49]. The radiated source has v = 0.5, and 
fo = 420 THz. The observer point is at x\ = 0.01, x 2 = 1.595, and the time is t = 2. The 
result of calculations is / = 417.82 THz, and the r = 3.1901. This point is indicated in Fig. 1(a) 
with / = 417.82 THz, and corresponds to v p = -0.31673, v g = 0.0084029, and x 2 < vt. Other 
solution was obtained for parameters x\ = 0, x 2 = 1.595 that results / = 428.9 THz and 
r = 3.1694 (see arrows in Fig. 1(a)). 



6 Conclusion 

In this paper the time-frequency integrals and the two-dimensional stationary phase method are 
applied to study the electromagnetic waves radiated by moving modulated sources in dispersive 
media. We show that such unified approach leads to explicit expressions for the field amplitudes 
and simple relations for the field eigenfrequencies and the retardation time that become the 
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Figure 2. (Color online.) Dependence of the Doppler shifted frequency u (u — 2wf) (a), and the 
time r (b) on the source frequency ujq for area where n(oj) < 0. See details in text. 

coupled variables. The main features of the technique are illustrated by examples of the moving 
source fields in the plasma and the Cherenkov radiation. In the paper it is emphasized that 
deeper comprehension and insight the wave effects in dispersive case already requires the explicit 
formulation of the dispersive material model. As the advanced application we have considered 
the Doppler frequency shift in a complex single-resonant dispersive metamaterial (Lorenz) model 
where in some frequency ranges the negativity of the real part of the refraction index can be 
reached. We have demonstrated that in a dispersive case the Doppler frequency shift acquires 
a nonlinear dependence on the modulating frequency of the radiated particle. The detailed 
frequency dependence of such a shift and spectral behavior of phase and group velocities (that 
have the opposite directions) are studied numerically. Such dependence in principle can be used 
also for reconstruction of unknown parameters of the dispersive medium. Other challenge is 
to enforce the developed approach to the Cerenkov radiation and its peculiar features in the 
dispersive and left-handed media. Such a problem will be considered in future study. 
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